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SUMMARY 


To understand fully the operation mechanism of a power semi- 
conductor device and in particular thyristor, it is important to know 
the minority carrier, electric field, potential and current density 
distributions within the device as a function of the spatial' variables . 
This report presents the formulation and closed form solution of the two- 
dimensional dependency of the minority carrier distribution within a 
thyristor for both the time independent and time dependent cases, 

Boundary conditions used to obtain the solution were chosen so as to 
conform to realistic devices being manufactured and to the physical 
properties of their surfaces. The general solutions were derived starting 
from the basic diffusion equations, and family of curves of the solutions 
for the static and dynamic cases of typical thyristors were obtained with 
the aid of a digital computer. These solutions were verified experimen- 
tally by the thermal mapping technique which was developed in the course 
of this work. This technique utilizes a 1.5y infrared beam to map the 
temperature in the device after the device is turned on for a certain 
length of time (static case) or switched on and off for a large number 
of times (dynamic case). The differential temperature in each point of 
the junction is measured by the attenuation of the infrared beam through 
the device whose temperature dependence was experimentally well determined 
Calculations were obtained by applying the analytical solutions to typical 
device geometries and making use of a CDC/3300 computer. The results and 
operating condition were plotted and compared to the results obtained with 
the experimental mapping method. The agreement between the analytical and 
experimental results were satisfactory. 


CHAPTER I 


INTRODUCTION ; 

A p-n-p-n device was invented by W. Shockley 1 in the early 50‘ s. 

The detaided device principles and the first working p-n-p-n devices 

• 2 

were reported by J.L. Moll et al in 1956. Later these works developed 

’ ' O 

into two categories: a) the p-n-p-n diode , or Shockley diode, is a 

classical switch which can change between a high-impedance "off" state 

and a low-impedance "on" state; b) the three-terminal p-n-p-n device^, or 

semiconductor-controlled rectifier (SCR), which has a gate terminal to con 

trol the switching, These devices are now developed to latch and control 

power systems because of their stable and low power dissipation in "on" 

and "off"- states. They are also used as dc choppers, inverters, frequency 
4 

converters , etc. , 

An extensive presentation of the basic theories and applications of 

O 

this device is available . The switching criterion was established and 
modified by others. T * * J.J. Ebers 5 established an approximate 
switching theory via a two-transistor model. From Figure 1 we find: 


! B1. = ^ V r A r C01 

(1) 

r C2 = a 2' l K + r C02 

(2) 

S&f" = I C2 

(3) 

= J A + h 

^ ( 4 ) 

we obtain: M Ml + ‘“l + 'c02 

A 1 - “1 - “2 

( 5 ) 



The gate current I G can trigger the switching, since when I Q is 
increased, I A increases, and the I vs. V characteristic is changed (Figure 

g 

2). When Ig flows into the Pg region, it will cause a non-uniform current 
density distribution which may over-heat, the device and cause abnormal opera- 
tion. This problem has been studied, * * ’ and many innovations of this 

device have been suggested. ^*13,14 An unc jerstanding of how the non-uniform 
current density is formed and changed during a series of switchings requires 
an understanding of the two-dimensional switching mechanism.^ It is the 
objective of this work to analyse the switching mechanism in two-dimensions 
and to illustrate numerically the process for typical operations. 

Experimentally, the lateral spreading of active-area (or turned-on area) 

9 1 fi 17 

was first tried by N. Mapham. W.H. Dodson et al , and H.J. Ruhr measured 

the spreading by measuring the currents of the probing islands (Figure 3). 

Multijunction semiconductor devices have rapidly become the main 
switching element for high efficiency compact power processing systems. 

Among these the thyristor is the most widely used because of its desirable 
means of control and latching characteristics and its superior ability to 
handle inductive loads. These devices must operate with the lowest possible 
dissipation in order to achieve high overall power processing efficiency 
under severe temperature conditions. Power switching devices and especially 
the thyristor are limited in their performance characteristics by the 
temperatures reached at the junctions. Nearly all significant operational 
parameters of a thyristor are affected in a major way by the temperature 
and thermal resistance of the junction. Also the life expectancy of the 
semiconductor power device is increased as the operating thermal stress is 


3 



reduced. When the device is switched on the current density dis- 
tribution within the junction varies greatly from point to point causing 
uneven temperature distributions. This could lead to local over-heating 
and become a source of hot spots and failure. It is thus essential that 
the thermal profile of the junction area of any power switching semiconductor 
device be known. It is not sufficient to measure the average junction tem- 
perature in order to rate the device, but it is required that the temperature 
at any point of the junction be measured. Such a novel method was developed 
on this project and is fully described in the Experimental Verification sec- 
tion of this report. In the 'work described in this report an optical system 
was set up to tes.t the temperature of a semiconductor device by measuring 
the infrared transmission which is temperature dependent The tech- 

nique will be extended to focus and scan an infrared beam on the device and 
thus derive the thermal profile within the device. The transmission varia- 
tions of the beam through the device produce thermal profiles which reveal 
the spreading of the active-area. 

In the following sections, the static analysis is discussed, Where the 
two dimensional time-independent minority carrier concentration will be 
solved for device configurations that have practical boundary conditions. 

The dynamic analysis is treated, where the time-dependent minority carrier 
concentration during the process of switching was derived for some type of 
device with very practical boundary conditions. Numerical examples calculated 
by digital computer are also given. The experimental work performed on this 
program is presented in some detail , including the description of the 
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ANODE 



ANODE 


(a) 





(b) 


I C = a *E*ho 
B = (l”Ct )Ig“ Icq 


Figure 1 

(a) Two transistor approximation of a three terminal p-n-p-n 

■ ; device..- vv-:; 

(b) Same as (a) using transistor notations 

(c) Current relationships in a p-n-p transistor 









laboratory set up, the techniques used in making experimental 
ments and in the analysis of the results. Also presented is 
of the results, their correllation, the conclusions and some 
tions for future work. 


measure- 
a discussion 
recommenda- 


CHAPTER II 


STATIC ANALYSIS 

In three-terminal semiconductor devices like transistors, thyristors 
etc., electric field, potential and current density are ndt simply func- 
tions of one spatial variable, since between emitter and collector, the 
third terminal always contributes lateral current disturbance. Therefore, 
it is necessary to use a two-dimensional analysis in order to investigate 
their operating mechanism in detail. For thyristors (Figure 4 and 5) the 
two-dimensional nature of the current flow is even more accentuated, since 
when the device is first "turned-on" by the application of a gate current, 
the load current density will start off slowly with a non-uniform distri- 
bution which will result in a partial turn-on of the device and could 

produce so called "hot-spots". This problem was particularly noticed in 

20 21 

the case of high-speed, high-power thyristors. * Two-dimensional 

approaches were used with limited success to derive the h-parameters for 

22 24 

bipolar transistors. * 

However, in order to determine the two-dimensional distribution of 
current, potential and electric field of a typical three-terminal semicon- 
ductor device, one must have a detailed knowledge about minority carrier 
distribution within the specified device and how this relates to its 
geometry. In this section an approach is employed that takes this into 
consideration and treats the general steady state two-dimensional soluti 
for the minority carrier distributions of semiconductor multi junction 
devices and the thyristor in particular. The novelty of this approach is 
that the analysis is applied to practical device geometries that are 
typically in use. 



Figures 


Figure 4. A typical thyristor geometry. 




Figure (5-a) A thyristor in forward blocking bias. 
Figure (5-b) Close up of p 2 - n^ regions 

Note, L " l 2 ~ £ 1 * anocle ’ 
k, cathode; G, gate 
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Formulation of the Problem 

For an n^ - P ? - n 2 transistor, a two-dimensional analysis should be 
applied to the P-region, and for a P-j - n^ - P 2 - n 2 thyristor, to the P 2 
region because that is where essentially the non-uniform power dissipation 
of the device occurs (see Figure 4). Starting with the continuity equation 
one can formulate the diffusion equation for the minority carrier concentra- 
tion n p 25 as 


D n v V X * y,t) + P n n p( x *y»t) v • E + n p (x,y,t)] 


= | ¥ n p (x,y, t ) + U n - G„ 


( 6 ) 


It is known, diffusion constant D is a non-linear, monotonically increasing 


function of n, 


26 


. p , but we may disregard this non-linearity if n p is less 
than the number of states in the conduction band Nc. 

To find a solution to equation (6) for the thyristor class of devices 
the following assumptions are made: 

(a) Geometry of the thyristor is as shown in Figure 5, and 
P 9 is the only region of concern. 

(b) The recombination rate of minority carriers U n can be 


represented by 
U • 


n P ' n Po 


(7) 


(c) Ionization generated by electric field E in P 2 - region can be 
neglected. The justification for this is that the resistivity 


of P 2 is about 1.2 ohm - cm, for a current density of 1 amp/cm^ 



the field E (< p j) is only 1.2 volt/cm. and the ionization 

C C 

generation rate of electrons is 3.8x10° exp(-l .75x10 /| E| ) 

?7 - 10 ^ -1 -3 . . 

|jj/q , or around the order of 10 sec -cm which is very, 

very small. On the other hand, for np - 10 cm" (which 0 is 

1 Amp/cm 2 ) and x n = 10" 6 sec, the recombination rate U n is around 

pi _i _•! . 

the order of 10 sec -cm , which makes the ionization genera- 
tion rate in equation (6) negligible. 

(d) No other carrier generations (e.g. surface generation) exist, 
thus G n = 0. 

(e) The oxide on the surface is considered as a perfect insulator. 
Namely, no current may flow through the oxide. 

(f) This analysis is limited to the low level injection case. 
Therefore, by charge neutrality, the Poisson's equation is 
now 

', 7.1* = 0. (8) 

£ 

and therefore the v . E term in equation (6) may be ignored. 

(g) The electric field E in the P 2 ~region is so small that the third 

term on the left of equation (6) can be ignored. The justifica- 

■ 

tion for high-power devices like thyristors, is that the field E 
is mainly in the x-direction (longitudinal). Considering the 
assumptions so far, equation (6) in x-direction can be rewritten 


for the static case: 

a 2 n^ 


’nO - * u n E 


8 n 


e_ = 


n P ” n P, 



ii 



The solution is 


n (x) = n p + C 1 
H o 


m-jX m 9 x 

+ C 2 e c 


where 


m 


• 2 ”I5 


i 


rt 

f.JL.’) + __ 

'2D ' D t 
n n n j 


(10) 


Thus the importance of E can be justified by comparing 
with 


4 D_ 


n n 


OD 


3 2 2 

For typical values of u n = 10 cmW-sec. D n = 25 cm / sec, 

T n = 10" 6 sec, E q = 10 V/ cm while the field due to ohmic drop, 

as mentioned in assumption (c), is 1.2 V/cm. The field due to 

doping profile in the P 2 -region is approximately (Dp/vipN^) 

(3N a / 3X), and is about 5 V/cm for a typical acceptor gradient 
18 -4 

of 2 x 10 cm . Therefore, as a first order solution, we 
may ignore the field completely in equation (6). As a result 
of the above assumptions, equation (6) can now be simplified 
for the steady state case as: 

n (x,y,t) - n po 


D 


p • V n p (x,y ,t) 


( 12 ) 


Equation (12) will be solved for the steady state condition for 
two distinct cases. The first is the steady state case with gate left open 
This corresponds to a forward blocking thyristor and forward conducting 
thyristor when triggered by gate pulse (current density in P 2 -region must 
meet assumption (g) ). This condition would also apply to a transistor 



with no base current. The second is the steady state case when the 
gate is in operation. This corresponds to the gate controlled thyristor 
(or typical transistor with base current). The transient case where the 
gate has just been applied will be treated separately in the following 
chapter. 

Substituting 

n s (x,y) = n p (x,y,o) (13) 


and defining 

N s (x,y) = n s (x,y) 



equation (12) may be rewritten as 

s v N s (x>y) = M x > y > 


(15) 


where, of course 


L = D t 
n n n 


(16) 



Boundary Condition and Solution 

Since for the forward blocking thyristor, the depletion region 
sweeps out all the electrons on the boundary we have rtp = o at x = o. 

(This is condition (A) listed below). Ohmic contracts force the boundary 
to maintain the carriers in equilibrium n pQ (condition (B) ). By assump- 
tion (e), current flow through the oxide boundary (i.e,, the concentration 
gradient normal to the boundary) is zero (condition (C) and (D) ). Con- 
servation of current should hold on the gates (condition (E) ) and junction 
J 3 (condition (F) ). When U 2 - ^) -* ■», we may have (see Figure (5-b) ), 


n s (W»2 ^ 




qV J3 / kT 



because in this case, P 2 - n 2 can be considered as a forward biased diode 
without "gates" (condition (6) ) . 

With reference to Figure (5-b) and equations (13) and (14) one may list 
the boundary conditions as follows: 


(A) N s (o,y) = n pQ 0 < y < L (18) • 

(B) N s (W,y) = 0 0 < y < £ 2 <y< L ( 19 ) 


■(G) 


CO) 


9N c 

1y 

3N. 


= 0 


ay 


y=o 


y< 


= 0 


0 < x < W 


0 < x < W 


( E > I nG -2q D 


r\ 


1 SN C 

W 


x=W 


dy 


( 20 ) 

( 21 ) 

( 22 ) 







w 




(F) 

■ . O 


f 9N c 

( 

3N 




q°„ H 

J 3X 1 

h 

dy + 

x=W J 

s 

3y 

y=^-] 

dx + 




w 



(23) 


(G) 

N s (Hi j) 

= % ' 

' n Po 

(&2 "T ^ ) 

• -> oo 

(24) 


9N C 

* 

( - 

v 3y 
y 

) dx 

'y=£ 2 


w 


NOTE: 


qO H ^ V J3 /kT 

l n ptt 3 ) * l~ n Po (e 


- 1) (2W - 2w + a 2 -A 1 ) 


(25) 


For forward conducting thyristor, junction J 2 is also forward biased 


(Figure 6), The minority carriers in proximity to the depletion region 
are no longer completely swept across J 2 - * n this case condition A, 
should be changed to: 


(A) N Q (0,y) = n pf - n p() 0<y<L 


where 


"Pf = a n Po e 


qV J2 /kT 


(27) 


Equation (17) should also be modified as 

qVjj/hKT 


n Pw = c n Po e 


(28) 


a.b.c, are the Interaction factors which account for the fact that each 


. . . oft 

junction is affected by the other two . Now, applying the separation of 


15 




Figure 6. A thyristor in forward conducting state; dashed rej^ons 
are depleted junctions. 
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29 

variables method, one sets 


N s (x,y) = X(x)Y(y) 


( 29 ) 


and equation (15) is then written, as 


32 v 32 y 

X — l + Y— 4| = XY 

3y 3x r 


■i 


(30) 


conditions (C) and (D) imply 

Y (y ) = cos j^- k - 0, 1, 2, ------ 


(31) 


By substituting equation (31) into equation (30), one obtains 


9 i 2 k 2 32 y 

l: 1 -x(\ k Y) + Y-i 

L 3x^ 


'n 


= XY 


(32) 


The particular solution for X is then of the form 


X(x) = si nh 



(33) 


We may satisfy conditions (A), (B), and (E) by superposing a few terms in 
the form of equation (33) for X(x) , so that N s (x,y) is -n^ 0 for x = 0, zero 
for x = W, and finite slope for x = W. A possible combination of terms 
giving a complete solution to equation (16) is as fol lows: 





(x - w)J (34) 

where U ■j(x-w) is a step function at x=w. B^, C k are unspecified 
coefficients for the condi tons mentioned above. 

As listed above, there are seven boundary conditions for this 
analysis, but only two (conditions (C) and (D) ) have been satisfied 
by the "cos part of equation (34). In this equation, there are only 
three sets of unknown coefficients, not enough to satisfy the other five 
conditions. One may superpose five terms in the form of equation (33) to 
satisfy these conditions, but then a complete solution with five infinite 


series will be too tedious to calculate numerically for any set of (x,y) . 
One may, however, simulate equation (34) to a certain function at a given 
x-plane, so that one brings in sufficient unknownsto satisfy these condi- 

m 

tions. One may assume a function for at x = w plane (see Figure 7) 
as .follows:.. 



0 < y < jl 

% < y < L -Jl 

L - sl < y < L 




Figure 7 


l 

1 



19 





n w» n o ,£ are t * ie new unknowns > in addition to A^, B k , C k , now one has 
the. exact number of unknowns for conditions (A), (B), (E), (F), (C), and 
equation (35) . 

3N 

The reason for assuming a function Q for -r— - v _ l( , instead of 

N s (w,y) is that it would bring faster convergence for the coefficients, 

and therefore reduce the number of terms to be retained in the infinite 

series. This is important to the numerical calculations to be performed 

The choice of "sin^e" for function Q was aimed towards a simpler mathema 

tical derivation for the coefficients as well as faster convergence for 

the infinite series. With the help of a digital computer, one may find 
4 

how close "sin 0 " is to an optimal choice. 

The determination of coefficients A^, and unknowns 

n“, n % j 1 should be related to the particular cases of interest, 
w o 

Case A I Gn = 0. 

(A) Forward blocking thyristor . 

By condition (A) and equation (34). 



n Po 

sinh(r— ) 

n 

0 


for k = 0 
for k f 0 





1 + 


for k = 0 
for k = odd 


(1 + 


W 


) for k even 
(37) 


where 


i i, 2 2 

r~ / 1 , k it 

VOt l—K + K— 

L L 

n 


From condition (E), 


C k = C k + A k cosh fa W 


COSh fa (W - w) 


(38) 


The details of these calculations are presented in Appendix A. 
By condition (F) ), 


I fjj * H ^ r^- (A„ cosh ~ + BJ 
np v 3 7 ^ n U n o Ln o' 


krrfi-i 

- 2 z A. sin — r— 
k even K 
■ kf 0 


Cosh v^w -r sinh 

Lv'cT 


kir 


(fa w)tanh va"(W - w) 
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Note, since 


and 


.fi n qV ,,/kT 

W J 3> ‘ 1“ 1)H(2W - 2w + ag - 


4 ) 

( 40 ) 



W + W 


(41) 


When, n po » p no> then I np (J 3 ) » I pn (J 3 ), and I pp (J 3 ) * 1 


np 


pn 


np 1 


K' 


By condition (G) and equation (34), 


n pw - n po E N s (w >2> 


= A sinh S- + s (-l)^ 2 A,, sinh VoT 
0 L n k i o k 


w 


k even 


(42) 


It follows from condition (B) that the minority carriers at the gate - 

metal contacts (assumed -to be ohmic) are at equilibrium. However, for 

a practical device the following two factors must be noted: 

1 6 —3 

(a) The p£ - base is usually lighly doped, (P 2 = 10 cm y 
typically), and the metal contacts can hardly be 
considered to be "ohmic" (true only for a P + - metal 
contact) . 

(b) The metal contact in this case is opened, and electrons can 
accumulate and build up a potential which tends to stop the 
sweep of electrons to the metal through the depletion region 


22 


of Pg - metal junction. This means that n s (lT,y) may become 

j, 

larger than equilibrium number n p{) , or N s (W,y) ^ n po for 
0 < y < ^1 , £ 2 < y < L, where 


N s (*r,o) 



cosh 


W - w 



E 

k o 
k even 



COSh /a(W - W) 

(43) 


By equations (36 - 39, 42, 43), we may solve for the unknowns 


V n o- 1 > V B k and c k- 


(B) Forward conducting thyristor 

Equations (37) to (43), obtained in part (A) remain valid but n pw as per 
equation (28) is replaced by n pw> From condition (A')» equation (36) now 
becomes 


' n po - n pf 
sinh r- 



for k = 0 
for k f 0 


(44) 


+ 

The force that sweeps electrons from semiconductor to metal by the field 
in depletion region can be balanced by the repulsion between the electrons 
in the metal and semiconductor. 

16-3 

For P« of 10 cm" , the maximum number of electrons on the metal (which 

£ C O 

is left open) is 2x10 cm" , then no electron can go farther than one 
third of the depletion length, c.f. The book by S.M. Sze (reference 25) 
Chapter 8 discusses the metal -semi conductor junction. 
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Case B 


w 0 


For a forward blocking thyristor, this represents the condition between 
gate triggering and turn on, or the condition the device is in after it 
has switched over. Equations (36) and (37) for and A^ in case A are 
still valid. By virtue of condition (B), 

sinh vcTW 

(45) 


c k = B k + \ 

sinh /a (W - w) 

From condition (E) arid equation (45) 


T Hji, A„ sinh r- 

Gn _ 2 n 

q D_ . W - w 

n sinh i 

n 

n 


kirJl- 


sin -- J-- ■ A^ sinh /a w 


+ 2 H l 


k even k 1 sinh ^ (W ' w) 


(46) 


and from condition (F) and equation (45) 

/ 


L - 1 % 


' H " °n < X L (A o cosh C + B o> 


kirS-! 

+ 2 E A. sin ~ r — 1 
k t o K L 


k even 


LVoT 

kir 


cosh, Sa w 


+ kfr sinh Sa w 


L7cT sinh voT (W - w) 


0 - 


cosh vcT (W - I 


( 47 ) 



The above equations are derived in Appendix B. By equation 
(36), (37), (42), (45), (46), and (47), one may solve for the 
unknowns n^, n', &, A fc , B k , and C k - 



Numerical Calculations 


Since the concentration of minority carriers in the forward 

12 -3 

blocking case could change from zero to 10 within 10 cm, the numeri- 
cal calculations may face a serious problem of accuracy. It was there- 
fore necessary to take particular care in the course of writing the 
computer program to minimize the possibility of errors. The calcu- 
lations were carried out on CDC-3300 Computer at Northeastern University 
Computation Center. The program was written in FORTRAN IV in double 
precision and the calculations were found to have an accuracy of at least 
twelve digits. The following examples are the numerical calculations 
corresponding to each case discussed above. 

Example 1 For a typical thyristor, in forward blocking state the 
geometry was as fol lows: 

H = 5 x 10" 2 cm, L = 5.1 x 10" 2 cm, w = 5 x 10" 3 cm, 

W = 1 x 1 O'" 2 cm, Ji-j = 3.5 x TO” 3 cm. 

-3 2 

The other parameters are, L p = 5 x 10 cm, D n = 25 cm /sec, and 

t = 10" 6 sec. Doping with P 2 = 10 16 cm" 3 , n 2 = 10 19 cm” 3 and 

4 -3 -7 

n po =2.2 x 10 cm . When I k = 4.4 x 10 amp, 1^ = 0, in forward 
blocking state one may solve equations (39), (42), and (43) simultaneously 
and obtain 

1.5118 X TO" 2 cm. 
n* = 6.5836652 x 10 13 cm” 4 ; 

: n* . = 2. 6683225 x 10 12 cm” 4 



Substituting these values In equations (36), (37), and (38) one 

obtains the coefficients A^, B^, C^, explicitly. Thus from equations 

(14) and (34), one could plot the minority carrier distribution 

n s (x,y) as shown in Figure 8. The derivative with respect to x, 
an (x,y) J nx (x>y) 

which is - ■ — . is shown in Figure 9. From Figure 8, one 

Vn 

can note that minority carriers distribute smoothly in P 2 -region. Figure 
9 shows how smooth the distribution of J nx is. However, one observes 
that in the region of w<x<W and 0<y<Jij , J px of y = 0 is larger than J nx 
of y = 2 .-| . This is because J n is mostly in the y-di recti on when elec- 
trons flow out of n 2 into P 2 region. 

• 

From this solution, n s (x,y), one calculates d n as 


J n = q D n^|x + yw ]n s (x ’ y) 


(48) 


One may estimate the electric field in x-direction as 


where 


~y -y 

E x = I px / [qup P 2 x(cross section area)] 


-y -y 


*Px “ 1 ’ *nx 


(49) 


(50) 


and 


I = H 1 J dy 
nx J nx y 


(51) 


Example 2 Same thyristor, but in forward conducting state without 
gate current. 1^ = 0.125 amp. Iq = 0, V 2 = 0.95 volt 
(polarity as defined in Figure 6), V 3 = 1.13 volts, 
a = 0.82, b = 1.5, c - 0.877. By equations (39) , (42), 
and (43), the parameters are calculated to be: 
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Figure 8 

The three-dimen- 
sional minority 
carrier distri- 
bution in P 2 " 
region, n p {x,y), 
in forv/ard block 
ing open gate 
case. To read 1 
the curves exa- 
ctly!, find the 
x- and y- axes 
first, then 
read on the 
vertical axes 
on each division 
of x-axis for 
the variables 
plotted. (Read 
the variables 
from y=o to 
y=jr ) • The 
curves i n 
y=L/2 to y=L 
are symertric 
to the curves 
shown here. 





V 
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n' = 1.35215768 x 10 19 cm" 4 
w 

= 2.01741901 x 10 17 cm" 4 

Again with these values, one obtains the coefficients An, Bn, Cn 

explicitly. The plot of n (x,y) is shown in Figure 10. The distri- 

3n (x,y) 3n (x,y) 
bution varies smoothly, - — - and — ^ - 

are plotted in Figure 11. From equation (48), one may find the 

distribution of J p in P 2 *-region. In the region of 0<y<^ , w<x<w, 

J n is flowing toward the gate when close to n 2 -region, but turned 

around and flows toward junction J 2 when close to the oxide (i.e., 

y = 0). This is due to the high current density at the P 2 - n 2 junction 

especially at a-j <y<£ 2 , for this forward conducting case. 






Figure 11a 

In forward conduc- 
ting open gate 


case, an p (x,y)/ax, 
the slope of 


n p (x,y) in x-direc- 


tion in p^-region 


(a). an p (x,y)/3y, 
the slope of 


n p (x,y) in y-di rec- 


tion (b). These 
two figures give 
the x and y com- 
ponents of 


v D „- 


an. 


ax 

(in cm' 4 ) 







CHAPTER III 


DYNAMIC ANALYSIS 


Introduction 

The basic principle of the switching of a thyristor can be well 

represented by equation (5) of Chapter I which is based on the assumption 

of uniform current density in the device. This is a one dimensional 

30 

approach which was amplified by E.S. Yang et al , and who pointed out that 

the switching is basically due to the current dependence of current gains 
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1 * 2 * J.F. Gibbons used a graphical method to analyse the I-V character- 
istics, but this is a time-independent analyses. The one-dimensional 
analysis valid for small area thyristors can be briefly states as follows: 
Consider the p-j-n^-pg portion of the device in Figure 12, the base 
current is 


! B1 = ! C2 
where I ^ is 
the ri 2 “P 2 "* n i 


the electron current collected by junction J 2 - 
portion of the device, the base current is 


(52) 

Consider 



+ I 


where 1^ 


is 


Cl 

the hole current collected by junction J 



(53) 

But,^ 


+ In a p-n-p transistor, the minority carrier distribution in base (of 
width Wg) can be approximated as, 

p(x) = p no (1-X/Wg) Exp(qV/kT) 

The total^charge of minority carriers is 

Q b = ^-fr- q P no Exp(qV/kT) 

where A is the area of cross section and 

I '* -AqD p (dp/dx) -- (2D p /W B 2 )Q B 




Figure 12. Model for turn-on analysis. 
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(54) 


Id - K,Q 


1 V 1 


J C2 " K 2 q 2 

From the equation of continuity, 


vol . of n 


/ 




9Q, Q1 

= X 


(-V. jp) dv = ^ + 


I 


So, 

! B1 


8 Q 1 Q 1 

3t T p 


likewise, 


. sq 2 , q 2 
B2 ' VT ^ 


(55) 


(56) 

(57) 


(58) 


(59) 


Substituting equations (54) and (59) into equation (53), one obtains 

(60) 


3Qo Q 

V¥i = *r + 7 


Substituting equations (55) and (58) into equation (52), one obtains 


¥2 ■ ^ + sr 


(61) 


Use Laplace transform to solve Q-j , Qg simultaneously in s^domain. 
then take inverse Laplace transform back to time domain. 
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Since 


I A (t) ■■ K 1 Q 1 (t) + K 2 Q 2 (t) (65) 

one may find the I A response as a function of Ig. 

1 2 

This approach was used by G.D. Bergman to calculate rise time and 
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delay time of switching and by R.L. Davies, et al . 

Instead of integrating the equation of continuity to establish the 
relationship between current (I) and charge (Q) as shown in equation (65), 
it is more rigorous to solve the minority carrier distribution from the 
diffusion equation, then calculate the electric quantities like current 
density J, field E and potential V. This approach assumes no lateral 
uniformity and applies equally well for the one- and two-dimensional 
analysis. This approach was used by A.A. Lebedev, A. I. Uvarov, and 
V.E. Chelnokov (one dimensional). 3 3>34,35,36 j^y solved the minority 
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carrier distributions in n-j and p 2 regions when "turned-on" by 
assuming no electric field, low level injection and that the number of 
minority carrier is negligible in "off" state. , They imposed the boundary 
conditions of the continuity of current at each junction and they obtained 
general solutions of minority carriers in the form of Laplace transform. 
They did not obtain the inverse Laplace transform and therefore did not 
obtain a general time dependent closed form solution and their time- 
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dependent solutions were presented only for some very particular cases. 

In the novel approach used in this work one solves the diffusion 
equation in n-j and p 2 regions separately for each region making use of 
the basic techniques like separation of variables, Fourier series and 
Laplace transform. These are general time-dependent two dimensional 
solutions applicable to practical thyristor device geometries. These 
solutions are then calculated numerically and plotted. With the computer 
programs given in Appendix C,D,E and F one may calculate the characteris- 
tics of other given thyristor types having different geometries and 
different doping characteristics. 



Formulation of the Problem 

As discussed in earlier chapter, the distribution of the minority 
carrier concentrations is a fundamental and descriptive approach to ob- 
tain electric quantities like J, E, V. This is true for both static and 
dynamic cases. For a typical thyristor one would start from the continuity 
equation and formulate the diffusion equations for the minority carrier 
P n in n^ region 


->- 1 2 -> i -* ->■ 

Dr> V p + u p V E-, *■ U Ev 

P F n M p H n 1 p 1 


p n" p no 


•1 

8p r 

9t" 


(V P„) 


( 66 ) 


where 


- 4 * 

v ' = 


9X 


+ y 


3y 


and for the minority carrier Op in p^-region 
D n ’ n p + “n "p " • h * y n h ‘ 0 "p> 


n„ - n„^ 9n„ 

_B Eo + _£ 

t at 

n 


( 67 ) 


where v = * W + y W 


E = x E x + y Ey 
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To find the solution to equations (66) and (67), one makes 
the following assumptions: 

(a) Geometry of the thyristor is as shown in Figure 13, 

Notice, the device is symmetrical along the 

y = L/2 line. 

(b) The recombination rate of minority carriers U and 

r 

U n are 

V p -A 

U n = V T " 

The equilibrium minority carrier concentrations of 
F no an< * n p 0 are negligible compared to P n and Op 
respectively for the dynamic case. 

. "-y - - * 

(c) Ionization generated by electric field E can be neglected. 
The justification is given under the section on Formulation 
of the Problem. 

(d) No other carrier generations exist, thus 
G p = 0, G n = 0 

(e) The oxide on the surface is considered as a perfect 
insulator. 

(f) This analysis is limited to the low level injection 
case. Therefore, the Poisson's equation is, 

^ -y 

v . E = p(x,y,t)/e * 0, by charge neutrality. 



-U -u w W 

J 1 J 2 


Figure 13 A Two-dimensional View of a Typical 
Thyristor. Slashed-line regions are 
depletion regions, potential polarities 
are defined as shown. 
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• ^ 

(g) Electric field E is not negligible in n-j-region, 

but is negligible in the P 2 -region. The justification 
for this assumption is also given in the section on 
Formulation of the Problem. In the n-| -region, one may 
compare the electric field (due to doping profile and 
current conduction) to the critical field to justify 
this assumption. As given by equation (11), the critical 
field is: 



2 

For typical values of = 350 cm /V- sec, 

2 -6 

Dp = 13 cm /sec, and Tp = 10 sec, the critical field is 

E Q = 20.6 V/ cm. The field due to doping profile in the 

n-j -region is (D n /y n N Q ) (aN^/ax) , and is estimated 

to be 2 V/cm and is negligible compared to the critical 

field. When the device is conducting a current of 
o 

1 Amp/cm (which is lower than typical forward conducting 
current) in the n-j -region of typical resistivity 40 ft-cm, 
the field is 40 V/cm, much larger than the critical electric 
field . Therefore , the el ectrf c f i el d i h n, -region becomes 
more and more important as the device is turned-on. 

(h) The problem is two-dimensional in P 2 -region, while in n^- 
region it is one-dimensional. The distinction is bused on 
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the geometry and whether there are more than two 
current sources or not. 

The boundary conditions can be stated as follows: 

»t) 

— Li 

y=o ' 9 y y=o 

y=L/2 y= L/2 

y=L y=L 




(A) n p (x,y,t) X=Q = p n (x,y 
t=o 


(C) n p (x,y,t) x= ^ = n pQ i 2 <y<l> 0<y<£ ] , 

r Jl = r J2 


'< E >' ■- ; rj2 + I fi- I j3 

Condition (A) refers to the fact that the number of minority 
carriers is zero at the depletion boundary, condition (B) refers to the 
insulation of oxide and symmetry of the device, conditon (C) refers to 
the ohmic contact at the gate contacts, conditions (D) and (E) refer to 
the conservation of current. Each current can be symbolically repre- 
sented as 






e qV-j/hkT 
( 68 ) 
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I J2 =LHM(Vj)q -Dp-^l x=u +„ P n (x--u) . 


I Sn n 

+ K(V 2 )qH 'o (D n ~5x x=o + “n n p (x=0 ) E 2x )dy 


iy an 
= <1 H ! 0l (°n “5^ 


n 17 x=w + y n V X=w) E 2x )dy 


w dn n 

+ 2 'w (D n1^ y=t 1 + u n n p^“ a l> E 2y> dx 
qV,/nkT 

t T ^ ^ 


+ l R3 e 




E-j , Ep are electric fields in n-| or regions respectively. 


E 1 ~ ^^1 ^ J1 ^ 


Vo o -1 

M(V 2 ) = 1 - ( - v - 2 ~ ) 3 

c. V n n o 


the multiplication factor, 


u ■ J 2 h *„ v A 


the depletion length of Jg. 


X R1 e 


qV^/nkT n.HL 


rA ^V 2 V e v 


qV,/2kT 
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the recombination current of 0^, and 

qV,/nkT n.HL / qV-/2kT 

J R3 e ’-T-Aw*** < 76) 

the recombination current of J 3 « 

Because the factor n (~2) is the exponent, the recombination 
current of J 3 is negligible in equations (68) and (70) when the deple- 
tion lengths of each junction are not much larger than the diffusion 
lengths. V gD2 for typical n-j (10 14 cm” 3 ) is 1150 volts, while the 
forward blocking voltage is no more than 500 volts, so equation (72) 
gives M(V 2 ) ~ 1. With these considerations along with boundary (A), (B) 
and assumption (g), we may simplify equations (68), (69) and (70) as: 


iji = L H q 


y p E lx P n^ x ”” u ^ " D p ax 


x=-U 


I J2 = L Hq 


-D. 


3P 


n 

ax 


x=-u 

+ qH °" V% 

' ! 

ax 


dy 

x=o 7 


(76) 


(77) 



r 

a n an^ 

w an 

1 

r 2 

dy + 2 / — & 
x-w J w ay 

<< 

it 

*=>■ 

Q. 

X 

■ JS| ax 

. . .. ■ 

L 


J 



We may obtain the minority carrier concentrations p (x,y»t) and 
n p (x,y,t) by solving equations (66) and (67) subject to the boundary 
conditions (A) through (E) and equations (76), (77) and (78). 
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Solutions for the Minority Carrier Concentration 

Equations (66) and (67) are two independent partial differential 

equations, but boundary conditions (D) and (E) joined with equations 

(76), (77), and (73) make p n (x,y,t) and np(x,y,t) mutually related. 

7 ft 9 

A. A. Lebedev et al. tried to solve them simultaneously * ’ ’ but it was 
too tedious to obtain the general solutions. We may, on the other hand, 
solve equations (66) and (67) independently, Ij^ , I J2 , and Ij^ are con- 
sidered as boundary values for ru- or P 2 ^region, the conservation of 
them (i.e. boundary conditions (D) and (E) will be the modification of 
these boundary values after certain interval . With the efficient 
digital computers available these days, the interval lengths can be 
chosen arbitrarily small, then the analytical and numerical solutions 
can be obtained with all the boundary conditions satisfied. 

(A) Solutions for n-j -region 

By considering assumptions (a) through (h), equation (66) can be 
rewritten as: 



In order to obtain the solution of this equation in closed form, it 
is found that to take the Laplace transform of equation (79) is a good 
approach. So, take the Laplace transform of both sides: 

" *1 “lx P ■ s, P - 

3 A 


where 


p(x,s) = ^f|p n (x»t) J = S 0 P n ^ x » t ) e " St dt 

a = : . u E-i/D 
pip 

s-| = (s + 1/Tp)/Dp 

p n (x,o) is the initial condition of P(x,s). It is the well-known 
static solution in the form of 

P n (x,o) ■ -p U sinh j—- / sinh 

P P 

= -b, sinh (81) 

P 

P n U = P n (x=-U) 
b-| s p a / sinh 

P 

The complete solution of equation (80) consists of a particular 

solution, p h (x,s) and complementary solution p r (x,s), 

P c 

P(x,s) = p p(x,s) + p c(x,s) (82) 

By substituting equation (81) into equation (80), the particular 
solution of this one-variable equation is 


bl c 

p p( x ’ s ) = 


(x+u)/L 


l+a,L -s,L 1 
Ip Ip 


-(x+u)/L 


s l L p +a l L p' ] 
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By substituting equations (82) and (83) into equation (80), we 
obtain the equation to solve for P c (x,s), 


V P c + a l £ P c ■ S 1 P c ' 84 > 

3 X 

P c (x,s) is obtainable by assuming it to be in the form of e mx . 
Solving form (two values), P (x,s) is obtained as 

V 

V a, 2 /4+s, (x+u) . -V77/4+Z (x+u) a, /v . llW9 

P c (x,s) - [ A 1 e 1 ... 1 + A 2 e 1 1 ] e 1 {x+u)/2 


Now, using the technique of inverse Laplace transform one obtains 
the time function p n (x,t). Using the transform 


,ht, _ , ht -s t 
i i a f e e 


Re {s> > h 


or, the inverse Laplace transform. 


Re {s> > h 


;,t) = 1 {P n (x,s)> 


b] (x+u+a-j D p t)/Lp 

v . n a,D t 

= - b ] sinh (^ + -!-E-) s - 
P L P ;V ; 
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Before obtaining the inverse Laplace transform, we would check 
the equation (85) first. By the "initial value theorem", 
lim s P\x,s) * p_(x,o) 

s>» n (87) 


Considering equations (82), (83), and (85), lim s P (x,s) 

S -> CD P 

and coefficient 0 f p c (x,s) has to be zero, since otherwise 

lim < p -A 2/ 4 + S 1 U+U) ■ - 

lim s e l 

s ■+■ 00 


P n (x,o) 


Note, x ^ -u<o in n-j -region. And also 


lim s e 

S -> OO 




74 + s-, (x+u) a 0 


Therefore, the value of is unrestricted by equation (87). As 
stated, boundary conditions (D) , (E) are now considered as boundary 
values. Condition (D) may then be restated as a description of grad- 
ient of p n (x,t). 


(D* ) 


3P n (x,t) 

ax 


X“-U ~ C 1 


which means Ip n (x=-u)=-HLqDpC^ * By a Laplace transformation, 

(D M ) - 8 - p (- x -» s ) = 

v ' ax x=-u s 


Condition (D") is not a rigorous condition, since in any time 
interval, the actual value of c-j varies as switching is started. There- 
fore, c, can be considered as a "final value" after a certain time 
interval , and 
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(D") i 1m s 


S + 0 


9£ 

9 X 


= Ci 


x=u 


Since 


9 X 


b-.L 

= JLL ( 

2 D. 1 


L_ ’ ) 

2 2 } 
l+a^p - s-j L a,L - l+s,L 


x=u - -p 1+a L 

1 h ~ n u p “Tp ~ ,T3 l t p 


+ A 1 ( ~~ + s 1 + a-j ^/4 ) 


(88) 


aj V 1 a " 

• • ' A 1 = C 1 / Sij- + —gr + ^~) 

/ P 


P c (x,s) = c 1 e 


-a. 


x+u 


V. 


/2 V 5 1 a l 

e' ST r? + T~ 
P P 


X+U 


s( !i + Vj_ + ^ 

sl 4 , 2 4 j 

P 


(89) 


Note, x+u = - 


x=u 


for x < -u 


-x+u = 


x+u 


tor x < -u 


By equation (58) and (59), and equation (10) of Bateman's Table 38 , 


1 { %TT~ } = e ' et [ e' ^ erfo^ /S- - j /ft ) 


+ e J ^®erfc (f /5~ + J^t )] 


( 90 ) 
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(91) 

From equations (86) and (91), the minority carrier concentration 
in the n^ -region is: 

P n (x,t) = P p (x,t) + P c (x,t) (92) 

To find the current one takes the x-derivative of p n (x,t) 

ap n (x,t) 3p p (x,t) 3p c (x,t) 

9X " " 3X + 3X 



cosh ( 


x+u 



a iV 

L 



P c (x,t) 



a-, (x+u)/2 


a l +Va l^ 4/L p 2 


-(x+uJZ^Zg 2 (x+u)Z T -Z 3 2 
e + e 
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/ 


^ + z, 


erfc(Z.j) e 


(x+u)Z, 


- erfc(Z 2 ) e 


■(x+ujZ, 


> (93) 


where Z 


I - ) 

1 s V 1/L p 2+a 1 2/ 4, Z 2 . -<*«)/ 2 + V t ''V a l 2 °p t/4 


Z 3 = -(x+u) 


/*vv- V 


t/ ! P +a i V /4j 


3P n (x,t) 

The solutions obtained for, p„(x»t) — r-— — , will be discussed 
later and numerically calculated for various locations and states of 
switching of the device. 

(B) Solutions for p 2 -region 

By considering assumptions (a) through (h), equation (67) can be 
rewritten as: 


J. 32 3n n 

V ,' — 7 + 7 ) n n = + 

n 9x dy* p T n 


( 


(94) 


Assuming that the variables are separable as discussed earlier, 
the function n p (x,y,t) can be separated as: 


n p (x,y,t) * 


? cos f|>(x»t) 
k=o L K 

k even 


(95) 


The choice of "cos " and "k" even is to satisfy boundary condition 
(B). 
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Putting equation (95) into equation (94), one obtains: 


2 ,.2ir2 an„ n„ 


3" k 


‘ a7'?' 'P 


( 96 ) 


or 


2 f 

a f k _ . f , . °'k 

— = af k + b 


8X 


at 


1 ,22 

where, a = — « + - — , , 

L ^ L- 
n 


b 5 


(97) 


Assuming: 

f k (x,t) = f k (x,o) + g k (x>t) (98) 

2 cos ^^-^(XsO) is the initial function of n (x,y,t), and is obtained 
as equation (34). g k (x,t) is the time-dependent part of n p (x,y,t). 

Introducing equation (98) into equation (97), one may separate 
time-independent terms and time-dependent terms. The time-independent 
terms have already been solved, and the unsolved part of the equation 
(the time-dependent term) is: 


a\ 


ax 


ag k 

? = a % + b W 


(99) 


We should solve equation (99) in the way that n p (x,y,t) meets 
conditions (A) through (E). For g k (x,t), condition (A) should be 


53 


modified as: 


(A* ) g k (o,o) = 0 

One may modify condition (A') further as 


(A M ) g k (x,o) = 0 


which means the time-dependent function g k (x,t) doesn't exist in a'ny 
place initially. 

Conditions (D) and (E) are now considered as boundary values, 
which are descriptions of gradients of n (x,y,t): 

■ • r 



where J n is the electron current density, J Qn is the electron current 
density at the gates. Gate current Iq is the current necessary to turn 
on the device, which can be turned on only when 1* is larger than a 
certain value. 

Since conditions (D ‘ ) and (E 1 ) involves time constants, we may 
rewrite conditions (D 1 ) and (E 1 ) for the boundary values of time-dependent 
function g k (x,t) in a time interval: 



8g k (x,t) 

sx 


x=o 
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(E") 


89 |.(x,t) 

ll \ In 


9X 


X=W 


" , "C”“ 

w 


c w =• 0 


Boundary values c , c- together with c, in Section (A) for 

Q W I 

n-j -region are fixed for certain time intervals, but will be modified 
at the end of each interval as one is checking with the conservation 
of currents (condition (D) and (E) ). 

Once again, the Laplace transform is found to be a good technique 
for solving equation (99) subject to conditions (A"), (C), (D"), and 
(E"). By taking the Laplace transform, equation (99) may be rewritten 
as: 


G k (x,s) = a s G k (x,s) 
9X 


( 100 ) 


where 

G k (x,s) s ^ jg k (x,t), 


a $ = a + bs 


The boundary conditions now becomes: 


(A") lim s G. (x,s) = 0 
s -> 00 K 


9Gi. 

(D") lim s 

S 0 


x=o 


c Q e'^“ <«>? 
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3G k 

lim S -rrr _ 
S V 0 ° x=w 



Condition (A") concerns the initial value of G k (x,s), while 
condition (E") concerns the initial value of 3G k (x,s)/9x at x=w 
which is proportional to the initial value of gate current Iq. Condi- 
tion (D") concerns the final value of 9G k (x,s)/9x at x=o for certain 
interval which is based on the same argument as in Section (A) for 
condition (D"). Parameter V is chosen so that G k (x,s) will carry it, 
and the choice of its value will ensure that condition (C) is satisfied. 
Expression "e "is chosen to ensure fast convergence when evaluating 

O 0 0 0 

the E series numerically (note, a = 1/L„ + k tt /L^) < 
k=o n 

k even 

The solution of equation (100) should be a combination of terms 
in the form of e ±>7 ^s^ w " x ^, e ±v ^s^" x ^, or e ±y/ ®s x . 

A form found appropriate for G k (x,s) is 
-»/a (w-x) 

G k (x,s) = A^ e - A 2 sinh vT x (101) 


which is mathematically easier than other forms for calculating the 
inverse Laplace transform of G k (x,s). A^ , A 2 , are coefficients to be 
determined. By defferentiating equation (101), one has 


dX 


r/a (oi-x) 

(A-| e - A 2 cosh /a" s x 


cosh /a w 


( 102 ) 
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where A 2 = ^2 CQS ^ w 


aG 


_k 

ax 




x=o 


= (A, e' s “ • A 2 / cosh ^s w > 


(103) 


aG 


_k 

ax 


x=w 


= /a s (A, e 


-v^ (u-w) - Ag ) 


(104) 


Condition (D") and equation (103), yields 


A, e"^ - Ag ' / cosh w = c 0 e'^ ' %Ji 


(105) 


while condition (E") and equation (104), yields 


A a' ’ c w e "^ /s/_a s 


Thus, 


c— e 
A, = ■-* 




— ~ 
s^ig cosh -/a w 


(106) 


and. 


A .2* 
M i 


c (7 


s/a 


s /a" cosh »^a w 
s 


(107) 


A^ and A 2 obtained are both positive, which ensures that G k (x,s) of 
equation (101) will always be positive. This is a necessary condition to 
have the number of electrons Hp(x,y,t) come out positive. 
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Now, 


G k (x,s) = ( 


s/a s/a s cosh /a w 


c- (“-X) 

) e 


— sinh /a x 

sVa_ cosh /a” w 


008 ) 


The inverse Laplace transform of this equation is based on 
the following inverse transforms: 


i, £ 


( \ 
-/a s (w-x) 


f Y 

-/s+a/b . v'b’ (us-x) 

e ... . \ 


/b /s+a/b 

( J 


i 
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Since a g = a + bs. From equation 5.6(6) of Bateman's Table , 


±' Id- 


/5" j /ift 


1 e ~a / 4t 


and 


£ 


G( s+to Q )/ 


= e 


'“ 0 t ^ 




'g(s)> 


the inverse transform becomes 


-at/b 9 

I - e e -b(o>-x) /4t 

^ /but 


0 09) 
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Also, 


i, - 


-1 A 




-^-1 J -Vs+a/b v'S (w-x) 


/ [(s+a/b)-a/b] 


= i [ e (“-x) ^ erfc (S^VI + Vp- ) 

» 


+ e- ( “‘ x) ^ erfc (^V“| . VI* ) 


(110) 


The last equality is shown in equation 5.6(6) of Bateman's Table. 
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Also, 


. £ y - 


nh Va" s x 


✓aT s cosh w 




(sinh /s+a/b ^ x) / /s+a/b cosh /s+a/b Sb w) 


_ p-at/b 6-1 ( — ) / bw 

' “ e 1 2U l*rb 


-at/b 

bW 




00 / ■. \T -bw^ (x/2w+r- t 

E '■'/ e 
^=-00 

(1U) 


The third equality is shown in equation 5.9(34) of Bateman's 

OO \ ' : ■' ' 

Table , 6|(v|t) is Theta function of first kind. The fourth equality 
is shown as the last equation on page 387 of the same Table. 

Note, equation (110) is the inverse Laplace transform good for the 
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first term of equation (108). Equations (109) and (111) are good 
for the second and the third terms of equation (108) respectively, if 
variable "s" is in the denominators of equations (109) and (111). By 
the property of: 

^' 1 *■ J'o s(x) dx 


and observing that the time functions of equations (109) and (111) 
are of the same form, the inverse transforms I-j , I 3 can all become 
directly applicable to equation (108) by: 

-aX-B/X 

I. = s* - e — dx a.6>0 

4 0 ft" 


2 s Jt e -aT 2 -B/T 2 


dt 


by x=T 2 



2V aj 

e erf ( /at + /g/t ) 


-2VaB 

+ e erf (/at /g/t )-2 si nh 2/a$ 


/a. /B>0 


012 ) 


The integration performed to give the last equality can be found 
as equation 7.4. 33 of HMF^. 

Therefore, the inverse Laplace transform of equation (108) can 
be obtained by I 2 and modified iy, I 3 when considering 1^, 



2& 
-/a(w-x) 


e (u-x)Va erfc ( y| + yat } 


+ e 
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2/E cosh /E w 
-/a(oi-x) 
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+ e erf(y?- - fj ) - 2 sinh //s"(u)-x) 


c w e 


-Wcu 
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2/E 


s (-1 ) r e 2^Sw|x/2w+r- 4-1 erf(/ /at + - 
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+ r 


2w 


+ e 


-2w//a |x/2w+r- j 


erf( yat. w ||_ + r4-| 4 ) 


- sinh 2 Va + r-^|w 


013 ) 


or, to make numerical calculations more efficient, equation (113) 
may be rewritten as: 


9 k (x,t) - 


1 


2/E 


/ a (w-x) 
e 


( c °' 


w 


_ erfc (Z ) 
cosh JE w ^ |J ' 


-^a (w-x) 

+ e (c + 


c w 


° cosh^ )( 2 - erfC( V> 
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where Z p ^ s /at/b + 0.5 (w-x) /b/t, 

Zp 2 s /at/b - 0.5 (ai-x) /57t, Z p3 = w x/2w-r~^| 

Note, when checking equation (113) or (114), we may find: 

Tim q^(x,t) = 0 
t o K 

1 i m 

t + « g k (x,t) - constant of time 

which means g k (x,t) is non-existing before switching and saturated 
after switching. 

Therefore, the minority carrier concentration in the P 2 ~region is: 

n p (x,y,t) = z cos [f k (x,0) + g k (x,t)] (115) 

k-o 

k even 

In order to calculate the current density of minority carriers 
in the p 2 -region, it is necessary to find the gradients of n p (x,y,t): 


62 


an p (x,y.t) 

ax 


- kirv 3f k (x » o) 

£ “ S T -It— 


k=o 

k even 




k=o 

k even 


( 115 ) 


The first term is calculated as equation (39), the second term 
can be calculated by: 


c- 


ag k (x,t) 1 

3x 2SK \ ' 0 cosh »/a w 


(c„ e' /?( “' X) C-^aTerfc (Z„,) e' Zp1 ] 
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0 cosh /a w 
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+c- e E (-1) (e p Ua erfc(Z 4 )-V ~ e p 


2/aZ 


f*= -oo 


-24Z n , r 
+ e p3 
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^(2-erfc(Z p5 » - e Zp5 


017) 


where Z p4 s /at/b + Zp3» // E7t, Z p5 = /at/b - Z p3 /b/t. 


an p (x,y,t) 

9y 


j {« s1n M 
k=o L L 

k even 


f k (x»o) + g k (x,t) 


(118) 


The solutions obtained, n (x,y,t), an (x,y,t)/ax, an (x,y f t)/ay 

n r r 

will be discussed later with numerical calculations for various locations 
or states of switching of the device, 
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Numerical Calculation of the Dynamic Case 
In the previous Section, equations (92) and (93) derived for 
the n-j -region and equations (115), (116), and (118) were derived for 
the p 2 -region. Now, using these equations, one may calculate numeri- 
cally the minority carrier distributions and current density with the 
aid of the following algorithm: (1) Establish the operating parameters 
of the device before triggering. (2) Determine the appropriate gate 
current I G to trigger the device. (3) Adjust the operating parameters 
to satisfy boundary conditions (D) and (E) for each time interval. 
Continue the above steps until the values obtained show a saturated 
condition, at that point the device stays "on' 1 . More specifically: 

(1) Parameters before triggering 

Set the device in the same operation conditions as in Example 1, 

_7 

for the static case, namely, 1^ = 4.4x10 Amp in the forward blocking 
case without gate current, which means Ij-j = I J2 = Ij 3 = 4. 4x1 O' -7 Amp. 
From the equality 

H L q D n p nri | qV-n/kT I n,. 

1)1* • no ^e J1 -ij coth (119) 

one obtains, 

p ■■ . p (e CV Jl /kT - 1) * Iji Lp / (HLqD coth p) (120) 

K nU K no v K K p 


Using the values of H, L, W, w» and parameters of L n> 

D n » x n from Example 1 of Section IV in Chapter II, (in addition, 

2 -3 -19 

knowing the parameters D =13 cm /sec.) L =1.14x10 cm, q=1.6x!0 

r r 

Coulomb). One finds from equation (120): 


P nU * 1.81xl0 23 x I 01 (120A) 

and, 

E, = Ij, / HLq p p p n - 8.8xl0 7 Ij, (121) 


From equation (119), = 0.69 V. 

In order to evaluate the switching characteristics of the device 
one has to connect this device in series with a load R , and apply DC 
power supply of voltage E v> The governing equation is: 


E v = IR c + 


\J1 


+ V J2 + 


J3 


(122) 


where I is the current that passes through the device. Assuming 
E y = 20 V, R c = 50 n, and = 0.69 V, V J3 = 0.6 V, (as calculated 
earlier), one obtains: 




V J3 = 18.71 V 


One may now evaluate the depletion length of u of n^ -region in 



junction J 2 : 


u * - 


J TO; = ' 


- x 10“ 4 cm 


(123) 


i-3 


u = 1.56x10 cm for this V J2 specified. 

Parameter c-j may be evaluated from p n (x,o) of equation (81) 


. m *n 
'1 sx 


x=-u 


p nU 

L p sinh(~^-T 


■1.17xl0 23 x I 


ji / 21 


nh 


U-u 


(124) 


-3 

Choosing U=3xl0 cm as a typical thickness of n-j-region, the only 
unknown parameters are then 1^ and V J2 . 

The parameters for the p 2 ~region before switching are specified 
in Example I of the static analysis- 
(2) Establishing a sufficient gate current 1^ 

A gate current 1^ is considered as sufficient to turn on the 
device if the time-dependent part of n_(x,y,t) is much larger than 

r 

7 fi 

the time-independent part, this occurs when t = lo to 10 sec. 
With I G specified, c w becomes, 


c w = 


3g k (x,t) 


3X 


x= , I* (125) 


As calculated by A.Y.C. Yu 4 ^, the minority carrier current I„ at 

bn 
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the gate contacts with semiconductor of 


p=10^ cm" 3 is 

I Gn =10 ' 6xI G 

Therefore, 

c = 1.15xl0 16 l r (125A) 

W u 


Numerically, a gate current is sufficient to turn on the device 
if c w (assigned by equation (125A) ensures the calculations of 
n p > 4x10^ cm" 3 and 9n p /9 x 2 5x10^ 3 at y=2.52xl0 -2 cm and t = 10"^ ~ 

ICf 6 sec. These numbers for n p and 9n p /3x are obtained from Figures 8 
and 9 respectively. 

(3) Parameter Modifications after Time Intervals 

When Iq is sufficiently large, I J2 and I J3 will be greater than 
Ij.j . After each time interval , one must then modify parameters to 
satisfy boundary conditions (D) and (E). From equations (76), (77) and 
(78), one may evaluate Ij-j , I J2 , I J3 . 
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18 |^280xE 1 
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-17 9 ^n 

I J2 = -2.856x10 
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|9 f (w,o) 8g n (w,i 
+ -fx— 


+ 2 
k=2 

k even 


L /_,„ kir . kia ] 

TciT (sin T - S1n ~ 


3f k (w,o) og k (w,t) 


z 

k=2 

k even 


ktT • 

IT S1 


ku A-i f ' v ' 

in lf k (x,o) + g k (x,t) 

'Ml 


dx l (128) 


When applying a sufficient Iq for a given interval, there are 
always three possible cases to adjust the parameters: 

(a) If Ij| is the largest among I J2 and I J3 , adjust c Q so that 

J J2 - T J1 and *J3 - ! J1 • 

(b) If Ij 2 is the largest among Ijy and I J3 , adjust Ij-j so that 
Ijl > Ij 2 > and adjust c Q so that I J3 I J2 . 

(c) If Ij 3 is the largest among 1^ and I J2 » adjust so that 
Ijl >_ Ij 3 > and adjust c Q so that Ij 2 1 Ij 3 - 

The "interval" lengths are chosen so that the differences between 
Ijl» Ijg and I J3 are below a certain fixed percentage. 

The numerical values of operating parameters can be found by 
using the algorithm mentioned above. (1) As 1^ = 1^ = 4.4X10" 7 
Amp. initially, then p n[J , , Vj-j , Vj 2 > u, c-j , are all known quantities. 
(2) With the computer program of Appendix C, one may find the sufficient 
gate current for switching to be, = 1x10"® Amp. and c- = 1.152x1 0" 1 0 
cm ^ . (3) With the computer program of Appendix D, one may calculate 
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the junction current at junction J 2 in p 2 -region as I J2 (n p ) = 

2x1 0~ 4 Amp, and C Q = 1x10^ 9 cm" 4 , when the gate current triggers 

the device for an interval of lo" 7 sec. With the computer program 

- 1 -7 

of Appendix D, one may obtain I J3 =1.1x10 J Amp after 10 sec. 

-2 

Also with the computer program one finds 1^ = 5.54x10 Amp, and 

_o 

the current at junction J 2 in P] -region r J2 (p n ) = 7.18x10" Amp after 

-7 -2 

lo" sec. Thus I J2 = I J2 (p n ) + Jj 2 (P n ) = 7 - 2xl °" Am P* Now ad J ust 

-1 22 

Ijl = 1x10 Amp, c Q = 1x10 for these programs and set interval to 
be 2xl0" 7 sec. and one calculates these values successively until the 
solutions for the condition of saturation are obtained. These solu- 
tions were plotted in Figure 14 for the n-|-region, and Figure 15 for 
the p 2 ~region. 


69 



■30 -18 


•28 


•24 


■22 


■18 


•16 

I 

J 

« « -u 


•12 


•10 


Figure 14 Minority Carrier Distribution when Triggering started 





Figure 14B x-derivative of p n (x,t) in the n-j -region 
when Triggering started (t = 10” 7 sec .) 
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Figure 15A Minority Carrier Distribution n p (x,y,t) when Triggering 
starts (t = 10~ 7 sec) (The way to read the curves is 
the same as that described for Figure 5.) 


2 



Figure 15B Minority Carrier Distribution n p (x,y,t) when 

the Device is Saturated in "ON" State. (t=0~ 6 sec.) 
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Figure 15D Minority Carrier Gradient Distribution np(x,y,t) 
(i.e. J nx ) When the Device is Saturated in "ON" 
State. (t-lO"® sec. ) 
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CHAPTER IV 


EXPERIMENTAL VERIFICATION 

The lateral spreading of the turned-on area for a thyristor 
(Figure 16), was reported by N. Mapham in 1962,^ who observed a time 
delay between the two gates when the device is turned-on at one gate 
and measured at the other. W.H. Dodson and R.L. Longini 16 subsequently 
measured the lateral time delay between probed islands using the con- 
figuration shown in Figure 17. The lateral voltage drop of forward 

23 

conducting thyristors was measured by E.S. Yang using a lumped 
simulating circuit as shown in Figure 18. In addition, as will be 
demonstrated further, since current density can be related to tem- 
perature generated in the device, an infrared microradiometer can 

also be used to map the thermal profile and give the current density 
42 

distribution. 

These approaches outlined above however give an integrated value 

41 1 ^ 

only, since the electrical measurements ’ * can not describe 

42 

differential values, and the microradiometer measurements read only 
the "smeared" temperature on the surface. Because the switching of 
such a device can be accurately described by the activities of minority 
carriers driven in and out of the p 2 -region, and because of the clear 
relationship of current density to temperature generated, it is possible 
to use the infrared technique to map the optical absorption and yield 
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Figure 16, 


aluminum 
contact 
gate 


A typical thyristor geometry. 
Areas shaded are ohmic con- 
tacts, the nonshaded surfaces 
are covered by oxide. Typical 
doping concentrations are: 


p ] = 10 

Po = 10 


19 -3 

cm . 

16 -3 

cm , 


nickel plated contact 
to cathod 
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✓ 
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cross section 

Figure 17. Device configuration. The island permitted the 
observation of the spread of the on-state. 


5 V 



Figure 18. A two dimensional lumped model for a p-n-p-n device. 








the thermal profile of the p 2 -region so that one may obtain the 

current density distribution at each point of the device junction. 

This in turn will be sufficient to describe the spreading of the turned 

15 

on area and the hot-spots generated due to successive switching. In 
this chapter the approach used in this project will be described, 
namely to employ a fine infra-red beam to map the thermal profile of 
semiconductor active devices. This is quite a novel technique. Infra- 
red micro radiometers were used in prior work, but only for surface 
measurements. This is the first time that an infra-red beam is used 
to scan and determine the differential temperature of an active device. 
From a knowledge of the temperature map, as will be shown, one can 
determine the current density distribution as well as other important 
parameters of the device. The development of the mapping technique 
itself used in this project is also novel and could be very important 
in determining hot spots, and device defects. The thermal profile of 
the p 2 region is based on the following two factors: (a) Infrared 
absorption rate of silicon depends on the temperature of silicon, 

(b) The major non-uniformity of current density is in the p 2 ~region. 
These factors will be discussed in the following paragraphs. 
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Theoretical Background of the Temperature Measurement Technique 

Near the absorption edge the optical absorption coefficient a can 
43 44 

be expressed as * 

a ~ (hv - E g )^ (129) 


for allowed direct transition, where h is Planck's constant and v is 

the frequency of incident light. The energy band gap of silicon E g is 

4-5 o o 

temperature dependent , between 300 to 600 k, with dE g /oT being about 

-2. 4x10"^ eV/°k. At 300°k, E g = 1.12 eV for silicon, therefore for an 

incident light of wavelenght X of 1.14 ym ( X = c/v = ch/E g ), we will 

observe a temperature-sensitive absorption coefficient. 

When the incident wavelength is less than 1.14 ym, absorption due 

to carriers being raised from the valence band to the conduction band 

becomes dominant; while when the incident wavelength is larger than 

1.14 ym, free carriers associated with impurities and defects contribute 

46 

the major absorption . In addition, the absorption coefficient of 
free carriers is 


a 


fc 


2gl_ 

A 2.,3 

4tt C n e. 
— o 


M n m n 




P P 




(130) 


Therefore carrier concentration dependence becomes explicit. In order 
to find the incident light to give the most temperature-sensitive absorp 
tion coefficient, it is necessary to find the spectral absorption coeffi 
eient of silicon wafers with various doping concentrations as a function 
of temperature. 
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For a small rectangular volume of silicon 
and resistivity p R , electric power dissipation 
volume is 

P = I 2 R * J 2 Pr AL 

Energy dissipated in time interval t is 

AE = P ( At ) = Pr AL J 2 (At) 

The equivalent heat is 

AQ = (aE)/ 4.1840 (in Joules) 

and the corresponding temperature change is 

AO 


where c y = 3 C^p 

n m = ~?8 ^ in moles ^ 
and Rq - 1.99 Cal/mole-deg 

Cpp = Dulong - Petit constant. The 

value varies with temperatute 

in °K^, between 300° to 500° 

with a temperature proportion 

—3 

of approximately 4.17x10 . 

m = mass = 2.42 for silicon 


of area A, length L 
of current I in this 

(131) 

(132) 

(133) 

(134) 


(gas constant) 
T 

k 9 

constant 


Therefore, in a time interval At, the temperature change in a 



small volume of silicon of resistivity p p by the power dissipation 
of current density J is 

AT- 1.66X10- 7 ( .pL jA t (135) 

l dp 

Note, the factor (p r /C dp ) is temperature-dependent, but its variation 
in the temperature range of 300° to 50G°K is negligible. 

Thermal conduction is a complicated problem because it is three- 
dimensional, involved with Fermi energy, mean free time, and differential 
thermoelectric power. The thermal conductivity is contributed by lattice, 
electronic or both lattice and electronic mechanism. It is related to 
doping concentrations'^ according to the relationship npy n Up/(ny n + 

pUp)^. In the temperature range of concern, namely from 300°K to 600°K, 

c 

the thermal conductivity for silicon with a doping level of 

n = cm" 3 varies from 1.6 watt/cm-°K to 0.52 Watt/cm-°K approximately 
51 

linearly . Therefore, the thermal conductivities can be characterized 
for this problem as: (1) Thermal conductivities for the p-j - or the 

n 2 -regions are one thousandth or less than that of or p 2 regions 
(2) The thermal conductivity dropped to 30% when the temperature 
raised from 300°K to 600°K for the n-j -region. For a block of silicon 
of thermal conductivity with temperature gradient [grad (T v ] the 
v heat flow through the boundary is 

Q th K th grad {T) (136) 
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Therefore the thermal conduction expressed by this factor is 
negligible in the or n 2 -regions with respect to that of n-j- or 
P 2 -regions. The thermal conduction is negligible in the n-j-region 
with respect to the p 2 -region as the temperature is raised from 300°K 
to 600°K. One may conclude, from the effect of thermal conduction in 
the p.|-, n-j-, or n 2 -regions, the thermal profiles formed in these 
regions are negligible compared to that in the p 2 -regions. 


Infrared Absorption Coefficients of Silicon as A Function of 

Temperature and Doping Concentrations 

A literature search for data on absorption coefficients as a 

function of temperature and doping concentrations for silicon revealed 

that the available information for this work is incomplete, especially 

for the doping concentrations of major concern for thyristors, i.e., 

18 -3 14 -3 

p - 10 cm , n ~ 10 cm . Therefore, it was necessary to conduct 
a series of experiments in order to determine the temperature and doping 
concentration dependence of absorption coefficients for various wave- 
lengths. 

A. Experiments in the wavelength range of 2.5 to 14 pm. 

A Double Beam Infrared Spectrophotometer can scan through the 
near infrared by modifying the reflecting grating, ranging 
from 2.5 pm to 14 ym. Basically this instrument (Figure 19) 
splits the beam emanating from a glowbar light source into 
two beams. They are then a.c. chopped, and one of them (which 
is 90° out of phase with the other) passes through the silicon 
wafer, the other bypasses it and then they are detected by 
the same detector, the ratio of transmission intensity to 
the incident intensity is available. A reflection grating 
rotating in the optical path in front of the detector allows 
the detection of a particular wavelength. A few things need to 
be noticed when operating with this system: (1) The water 

vapor band around 6 ym disturbs the readout. This can be 
corrected by flushing nitrogen into this system to clear the 
water vapor, (2) The beams are around the size of 5 mm, 
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BEAM 

SAMPLE COMBINING 

SOURCE AREA AREA AREA 



MONOCHROMATOR 


DOUBLE BEAM FAR INFRARED SPECTROPHOTOMETER 


Figure 19. Optical Diagram 




and can not be focused without disturbing the whole 
alignment. So, it is not suitable for use in mapping the 
thermal profile. 

A sample holder was made to hold one-inch silicon wafers 
(Figure 20). A thermocouple was placed in contact with the 
wafer so as to have good control of the temperature of the 
wafer. Quartz was used for the windows because it has 
negligible optical absorption for the temperature and wave- 
length range of concern. The wafers used were cleaned and 

polished 0.02 cm thick silicon single crystal wafers with 

16 3 18 3 

doping concentrations of p = 10 cm , p = 6x10 cm , 

14 -3 18 -3 

n = 5x10 cm , n = 8x10 cm . Temperature variations 
were sampled as 25°C, 90°C, 145°C, and 220°C. The trans- 
mitted intensity changes relative to the intensity at room 
temperature were measured. This set of data shows that the 
intensity change varies with wavelength and with doping con- 
centrations. For p-type silicon, major changes are around 
3 um and 12.5 pm, while for n-type silicon, major changes are 
around 4 pm only. None of these intensity changes give more 
than 10% range for all four types of silicon, which implies 
that this data can not be used to resolve temperature profiles 
satisfactorily when the wafer temperature changes from 25°C 
to 250°C. 

B. Experiments with wavelengths 1.15 pm and 2.0 ym. 

An optical system as in Figure 21 was set up to measure 
optical absorption coefficients for monochromatic light. 
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i Wafer 


Thermal Couole 


Power Supply 


Heating Coll 


Figure 20 


Sample Holder with Temperature Control and 
Probinr Thermal Couple 
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Figure 21. Optical system to calibrate silicon optical 

absorption coefficients versus silicon doping 
concentration and temperature, 
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Tungsten light was used as a source, the beam was chopped 
at a frequency of 80 cps and samples were held in the tem- 
perature controlled holder. The 2.0 ym filter used was manu- 
factured by Optics Technology Inc., Flushing, N.Y., the wave- 
length tolerance of this filter, set 20 Infrared Narrow Band 
Filter, was ±0.05 ym. The photovoltaic indium arsenide 

infrared detector, when cooled by liquid nitrogen to 196°K, had 

10 JL. 

a peak detectivity D* of 3.7x10 cm-Hzvwatt at 3.3 ym, 

the output impedance was 3000 ft and a response time of 1 ysec. 
This system was checked by varying the intensity of the source 
from maximum to 5Z of maximum and the oscilloscope does show 
a corresponding amplitude variation of the signal with errors 
below 6%. 

Data were taken on sampling wafers of various doping con- 
centrations and on temperature settings of 25°C to 250°C. 
Decrease of transmitted intensities versus temperature were 
measured and plotted on Figure 22, for the incident light of 
1 .15 ym. 

From the relationship 

I = I exp (-ad) (137) 

one may find a. I and I are the intensity of light before 
and after passing through the sample of thickness d, and absorp- 
tion coefficient a. The a's calculated are plotted on Figure 23 
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Figure 22 A. Transmission change of infrared light (\-),15um) 
as a function of temperature for silicon wafers. 
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Transmission change of infrared light (x=l 
as a function of temperature for silicon wafers. 
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Figure 22C. Transmission change of infrared light (,\=U5ym) 
as a function of temperature for silicon wafers. 
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Figure 23C. Absorption coefficient change of infrared light (\=1 . 15pm) 
as a function of temperature for silicon wafers. 
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Figure 23D. Absorption coefficient change of infrared light (A=1.15ym) 
as a function of temperature for silicon wafers. 




In the same way, incident light of 2.0 urn, the intensity 
changes were measured and a's were calculated. 

The characteristics of Figure 23 plotted for a's are: 

(1) the higher the doping concentrations, the larger the 
a's. (2) The higher the temperature, the higher the a's. 
These findings are consistent with the theories given earlier. 
In fact, these data are also consistent with the experimental 
data of references 45, 46, 19. By comparing the set of data 
of a's for x=l .15 ym to the set of a's for A=2.Q ym, it 
appears that a's measured by the light with a wavelength of 
1.15 ym can resolve temperature better than a's measured at 
other infrared wavelengths. The optical absorption coeffi- 
cient has the highest temperature sensitivity when the 
Incident light has an energy value around the band gap energy 



The Mapping of Thermal Profiles 

Using the data obtained in an earlier section, one may find 
the temperature of the wafer of known doping concentration and thick- 
ness from the a - T curves obtained. Therefore, by scanning light 
through a device and measuring absorption coefficients continually, 
one may draw the thermal profile of it. The problem involved with 

this experiment are: (1) The cross section area of thyristors are 

-2 2 

usually around 10 cm or less, the incident beam spot diameter 
should be around 2.54 ym, so as to have enough number of discriminat- 
ing measurements of absorption coefficients. (2) The device should 
not be covered with metal or thick oxide over these regions where light 
is transmitted. (3) The scanning of the device should be well 
arranged so as to enable each measured point to be identified. These 
problems are the major objectives in setting up the optical system as 
described below: 

A. The Components of Optical Scanning System 

Everything in the system should be solid, well-aligned 

and tightly leveled on an optical bench (see Figure 24): 

Each component is specified as follows: 

(1) Light source: A He-Ne laser was used as Tight source, 

its cavity length is 43 cm and designed to emit light 
o 

of 6328 A in TEMoo mode with power of 2.0 mW. With 
confocal mirrors it emits light of wavelength 1.15 urn 
with 0.3 mW power. Amplitude was found stable with 
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Figure 24 Optical System to Map the Thermal Profile 
Of a Switching Device 




2 

ripple within 1%, beam diameter is (1/e point) 

1.5 mm. and beam divergence is 1.5 inr, and linearly 
polarized. 

(2) Chopper: It chops laser beam to square waves of 

frequency of 80 cps. 

(3) Lens A: An epiplan 40/0.85 lens was used. 

(4) Device holder: An optical bench holder with micrometers 

for adjusting the device location horizontally and 

-3 

vertically with a precision of 10 cm, was found adequate 
to localize the focused spot. 

(5) Lens B: The focal length is 9 cm which refocuses the 

diverging beam on the photodetector after the beam is 
reflected by the front surface mirror. 

(6) Photo detector: The same InAs photovotaic detector 
was used as was used previously. 

Devices: 

In order to verify experimentally the actual temperature 
profile within the junction, it is necessary to test typical 
Commercially available thyristors. It is however very 
difficult to know the doping concentrations, junction depths, 
metalization layers, etc. of a commercial product. The set 
of measured figures one can obtain is that shown in Figure 25 
(GE Series C30 which is corner-gated). When measuring along 
a-b line or c-d line of Figure 25-B, the activities along a-b 
line are too far away to influence the c-d line, and vice 





Figure 25-B Geometry of Gate-Cathode Metal Contacts 
of Thyristor under test. Solid-lined 
regions are the new contacts after 
modification. The anode on the other side 
is modified to the same solid-lined 
geometry also. 






versa, so that the measurements of a corner-gated 

thyristor should have the same characteristic as an 

interdigitated thyristor. This Series C30 thyristor has 

4 

a power switching capability of 10 Watts, maximum forward 
blocking voltage of 400 volts, RMS forward current (on-state) 
25 Amperes, gates power dissipation of 0.5 Watts (average), 
or 5.0 Watts (maximum). The metal contacts are made of three 
layers: chrome (10 3 X), nickel (3 ~ 4x1 0 3 A), silver 
(1.5x10^ H). Since the metal contact on the anode is a 
complete 0.38 cm x 0.38 cm square, it blocks the infrared 
beam. This was solved by modifing the metal contacts of 
this device on a solid holder and attach leads to it is a 
very difficult task, because the device has no direct sup- 
port and numerous attempts to bond wires with ultrasonic 
wire bonding were fruitless. An alternative approach was to 
take a piece of printed circuit board (2 inches in diameter) 
cut a hole in the middle and glue a chip of devices onto it. 
The copper surface of PC board was cut into three areas and 
spring contacts were attached to the anode, cathode and 
gate of the device from each area individually as shown 
in Figure 26. The modified devices were tested with a Tek- 
tronix Type 575 Curve Tracer and the I-V characteristics 
obtained are shown in Figure 27. 




Testing Thyristor mounted on Printed Circuit 


Figure 26 
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Figure 27. I-V characteristics of modified 
sample C30. 




C. The Set-up and The Measurements 

To set up a system as shown in Figure 24 with components 
specified in A. is fairly straight forward, but to make sure 
that the beam size is small enough, a pinhole was used in the 
optical path of this system. The pinhole chosen was 25 ym in 
diameter. The position of the pinhole is set to obtain a 
maximum intensity detected (the CRT reads 0.07 volts for 
peak to peak magnitude of the square waves, and when measured 
without the pinhole, the CRT reads 0.14 volts). The spot size 
focused by the microscope objective as described in IV A (3), 
was about 0.00375 cm in diameter. 

A photograph of the actual experimental system is shown in 
Figure 28. The Curve Tracer monitored the device characteris- 
tics during measuring, and it provides switching drive at a 
frequency of 120 cps. 

When measuring along the ab line or the cd line of each 
device, the thermal activities may be mapped between 
cathodes and gates. There are two kinds of thermal pro- 
files to generate: (A) Thermal profile of the device in 

on-state, (static case). (B) Thermal profile of the device 
due to successive switching. Since the temperature change 
due to each switching is accumulated, this thermal profile 
may reveal the change due to each switching (dynamic case). 
Dynamic measurements were conducted first. 
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Figure 28 The Actual Experimental System 
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(A) Dynamic Measurement There are two ways of measuring: 

(1) Scan the beam through the device (in fact, one moves 
the device, so that the fixed light may "scan" through the 
device) before the device is energized, and take several 
measurements of transmission during that time, then switch 
the device on and off a large number of times, scan the 
beam through again, and measure the transmission at the 
same spots as the first time. This should reveal the 
difference of transmission due to switching. But this 
approach needs to synchronize the system and be fast 
enough to avoid thermal runaway in the second measuring. 

(2) Fix the beam at a spot on the device, measure before 
or after the switchings, then cool the device to room tem- 
perature and move to the next spot. Since it is difficult 
to take care of the thermal runaway in the first approach, 
the second approach was tried in this research. When doing 
the measurements, the scanning lines (i.e., ab or cd lines) 
were 0.01 cm away from the edges. The spots for measurements 
were not chosen equidistant. Between two measurements at each 
spot, the devices were switched on and off for five minutes (a- 
bout 3.6x10^ switchings). After the second measurement, the 
device was cooled down by an air blower until the intensity 
went back to that of the first measurement, this serving as a 
consistency check. The results of these measurements are 
tabulated in Table 1 through 5. Column "y" identifies 
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distances between measuring spots and the metal contact 
of the cathode. Column "V-|" is the intensity before switch- 
ing and "V Q " is the Intensity measured without the device. 

"a" is the temporal absorption coefficient immediately after 
switchings, and "T" is the corresponding temperature. 

(B) Static Measurement The set-up is the same as that shown 
in Figure 24, only the curve tracer is replaced by a DC power 
supply. The positive output of the power supply is connected 
to the anode of the device. It is also connected to the gate 
through a load register Rg. The negative output is connected 
to the cathode through another load resistor R c - When the 
device is turned on for a moment, a non-uniform thermal dis- 
tribution is generated. Here again the same two methods of 
measuring the thermal profile which were described for the 
dynamic measurements are available, and the second way was 
chosen. For these measurements, Rg was chosen as 200 a, 
and Rg was chosen as 10 ft. The power supply was set to 8 

volts, when in forward blocking state and the cathode current 
-4 

1^ was around 5x10 Amp. When triggered by gate pulse, the 
device was turned on, 1^ jumped to 0.125 Amp, and the device 
voltage drop was 0.2 volts. Optical measurements were 
taken before switching and thirty (30) seconds after they 
have been switched on. Then the device was turned off by 
turni ng the power supply off. 
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Absorption Coefficients and Ter.rerature Values after Thyristor wn-s 
switched fcr ^.6x10** tir*ies. V ■* 0.14 volts. Thyristor 5 



The device was cooled and the beam shifted to another 
point thus preparing it for another measurement. The re- 
sults of these measurements are tabulated in Table 6 
through 10. The caluclated absorption coefficient and 
corresponding temperature T are also shown in these tables. 
D. Data Analysis 

(A) Dynamic Measurements From the data for each of the 
points on the device, the absorption coefficients at a 
given temperature can be calculated and interpreted. Fol- 
lowing is a typical calculation: From Table 1, 

when y = 42.5, V Q = 0.14, V ] = 0.045, 

V 2 = 0.042. From Figure 23A, a Q = 47 cm" 1 for 
16-3 

p = 10 cm at room temperature. Since the measure- 
ments were done on the gate contact region, the thickness 
of p 2 is 45 pm (instead of 30 pm for the typical P 2 re- 
gion). 

The transmission ratio is: 

I _ - ad _ -47x0.0045 _ n P1 

t - — e o = e =0.81 

x o 


But V-,/V 0 = 0.322 


Therefore the transmission ratio due to p| and ny is 


0.322 


0.395 


After switching on and off for five minutes, V 2 is ob- 
tained, and V 2 /V Q = 0.300. the actual transmission ratio 
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due to Pg» when considering 39 . 5% of transmission is 
due to p-j-n-j regions, is 


0.300 _ 

055" 


0.76 


This 76% is due to the new a and the same d(45 m) of 

P 2 -region, i.e. 

-ad -0.0045a n 
e = e =0.76 

. ’ • a = 60.5 cm 1 

From Figure 23A, one finds that the corresponding 
temperature is 178°C. 

The same calculations go on to complete columns of 
Vj/V', Vg/V^a, T of Tables 1 to 5. The average values 
of a and T were plotted in Figures 29 and 30. 

Note, (1) The transmission factor 0.395 due to p^ and 
n-j is consistent with the data measured. From 
Figures 23B, for p-^ -region, a = 172 cnv" 1 at room tempera- 
ture and a = 85.6 cm" 1 for n-j -region at room temperature 
(Figure 23C). From the thickness of p^n-] given in 
Figure 25-A, the transmission is 

e x e = 0.389 

which is only -1.5% off 0.395. (2) From Figures 23B and 

23C, it is apparent that the p-j and n| regions are much less 
temperature sensitive than the p 2 -regions. Suppose, after 
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many switchings, there is a hot spot of 240°C in p 2 , 
the neighboring point in n 1 is at 160°C, and p 1 is at 
80°C, the transmission is: 

-182x0.045 v -87.2x0.0025 n „ Q , 
e x e = 0.386 

which is only -0.8% off 0.389. However, the transmission 
ratio in P 2 goes from 

e -47x0.0045 0 81 

to 

e -67x0.0045 . 0< „ 

which is -8.65% off 0.81. Noting the difference of 
thermal conductivity of each region as is mentioned 
earlier, it follows that the thermal profiles obtained 
in Figures 29 and 30 are mainly generated by the p 2 
region. 

(B) Static Measurements The analysis is essentially 
the same as that shown for Dynamic Measurements, the 
values of a and T are calculated and tabulated in 
Table 6 through 10. The average values of a and T of 
these five measurements are plotted in Figures 31 and 
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Variation of the temperature value between cathodes 
and gates. (Dynamic Case). 
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Figure 31. Variation of the absorption coefficient value between 
gates and cathodes. (Static Case). 
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Figure 32. Variation of the temperature value between cathodes. 
(Static Case). 
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CHAPTER V 


RESULTS AND CONCLUSIONS 


The static measurements shown in Chapter IV Tables 6 through 10 
and Figures 31 and 32 may be compared with the theoretical calculations 
of Chapter II. From the current density-temperature correspondence 
given earlier in equation (135) namely, 

A T * 1 ,66xl0” 7 ( -!X ) J 2 At 

l DP 

one sees that current density J is proportional to //aT. The current 
density of minority carriers in the p 2 -region as calculated in equation 
(35) is of the form 


( + (K, “ n 0 ) si n 4 |f- 


9 N 

J np 2 = qD n Tx 


* 


w 


0<y<Jt 
«.<y<L-£ 

+ (n w - nj$in 4 L-a<y<L 


w 


2T 


and is plotted in Figure 7. 

One could compare these two results by fitting the expressions 
of "Sin 4 to the plot of a/aT\ This is done in Figure 33, by choosing 

the value of 'V in the expression of Sin 4 to be 2x1 0" 2 cm. 
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The experimental result ofa/TTgives an overall thermal profile of 
the device, while the analytical calculation show a localized current 

density distribution J, e.g. equation (35) of Chapter II is giving 
aN 4 

at x=w. The assumption of Sin H in Chapter IT fits the experi- 
mental results satisfactorily when V is well-chosen. At the same 

I I 

time one can obtain n w and n Q to specify the "central" and "outer" 
current density. The thermal profile measured and shown in Figure 32 
of Chapter IV, orVAT shown in Figure 33, is matching with the boundary 
condition 


- 3 

In the range of y=(50' , 80)xl0 cm the curve does not quite stay level 

as suggested for the central region of Figure 7. This may be due to the 

crystal defect and non-uniform doping depth that cause a non-uniform 

52 

resistivity for the device in the "on" state . Also, when modifying 
the contacts, the photoresist processing is not smooth enough to give 
well-defined metal contacts, and this could account for the level 
variations. 

The dynamic measurements shown in Table 1 through 5 and Figures 
29 and 30 may be compared with the theoretical calculations. By the 
linear correspondence between VTT and J in equation (135), the measured 
AT of Figure 30 should be plotted asWAT in Figure 34, and typical 
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current density calculation in Figure 15D plotted on this same graph 

I - - ' ; _ 

and the two compared. When this is done one can conclude that: 

1. The correspondence between analytical and experimental 
results is quite satisfactory. 

2. In the current density plot in Figure 15D, the curve that 
fits the experimental result best is the curve that represents 
the dominant current density distribution in the interval of 
triggering. 

3. The mathematical methods developed are useful in obtaining 
solutions sought in this work. 

4. The current density distribution assumed for the static 
analysis was valid. 

5. The infrared thermal mapping technique developed and used in 
this work produced thermal profiles, which verify the static 
and dynamic analysis. 

This technique could also be very well used to identify "hot spots" 
of thyristors or other multi junction devices. The analytical results 
verified serve as descriptive functions for minority carrier distribution 
at all points at each instant of time. The method used to calculate 


these functions numerically is a valid method of describing the switching 


and will help understand the operation of these devices. The analysis 


presented in this report can be further extended to some special case, 
like high level injection. It can also serve the device designers when- 


ever numerical calculations are done with various device parameters and 
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operation conditions. The mathematical techniques used in the 
analysis is not limited to thyristors, but are also applicable to other 
three-terminal devices like bipolar transistors. 

The results derived from this work have made a significant contri- 
bution to the methods of analysis of thyristor (and other multi junction 
type devices) and have given good verification of these methods by ex- 
perimental means. An equally significant result of this work is the 


development of the experimental method of testing these devices for 


design optimization hot spots, fabrication defects and improper device 


operating setup by examining the temperature map within the device. The 
experimental setup can be further modified so as to build a compact 


system to map thermal profiles of semiconductor devices. These new 


extensions of this work are well worth pursuing. 
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APPENDIX A 


STATIC CALCULATIONS OF A R AND B k WHEN I = 0 


By condition (A) and equation (34), we have 


n_ * £ cos B. sinh /a w 

P° k-o L k 


where a means 




1 + kV 

L 2 L 2 
n 


Integrating both sides with respect to y from 0 to L, we obtain 

( 

n 


JEO. 


for k = 0 


\ = ( s1nh c; 

0 for k f 0 

Taking the drivative of equation (34) one obtains: 


V 


5N e 

W 


= 2 COS /a [A. COSh /fit W + B. ] 

x=w k=o . K K 


Equating this to equation (35), one can find the values of A^. 


u / L c ° sh r + B o> d * ' f 

n*' o n vo 


dy 
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(A-T) 

(A-2) 
For k 


For k/0, 


A o * | K< L ■ I ^ * 4 V l" 


- B 0 !>/ cosh r 


(A-3) 


A k cosh /a w 


A 

I ... k TT 

J cos “T 


X- ( £ cos )dy 

Mo L 


-L-A 


z 1 


n w i cos dy + J + J n 0 cos dy 


o L-a 


+ 1 K - %> sln l?F cos x*- 


dy 


+ / (% * n^) sin 4 ir(l - y) 
L-Jl w 0 2A 


2 F - cos d y 


(A-4) 


and 


A,. = 


4(n; - n-) sinfe 
ntr r/a cosh n/a W 


- 1 + 


k£ 


L(2 + j--) 


1 + -Sr 0 + 


1 + 


ka 


i + -,,-L - + 1 


2 |i . 4 2 M. . 2 


k's are even and non -zero integers. 
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APPENDIX B 


STATIC CALCULATION OF A K , B K AND WHEN I/O 


When I R / 0, the gate is closed, no charge can accumulate on the metal 
contact which will block electrons from sweeping out, thus condition B, 
is now valid and 


N (W,y) = £ cos [A. sinh (/a W) + B, sinh (A" (W - w)) - C, 

s k=o L K K * 


>/a (w - w)) 


= 0 


(B-l ) 


In order to obtain A^, and C^ independent of y, we may set the contents 
of the square bracket to zero, namely. 


Ci. = B,. + A 


sinh^a W 


k k k sinh[/a(W - w)] 


(B-2) 


By condition (E) and (B-2), 


W ■. 2 Ht 1 A 0 s1nh(W/L n ) 
q D n l n sinh[(W - w)/L n ] 


+ 2H I 
k / 0 
k,even 


/at sin(k7rA^/L)A k sinh(yaw) 
for sinhE/cT (W - w)] 


(B-3) 
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APPENDIX C 

PROCESS TO MODIFY METAL CONTACTS OF METALIZED THYRISTOR 

' I ; 

The process to modify the metal contacts is as follows: 

(1) Strip metal contacts: Put a wafer of finished but uncut 

thyristors into ferric nitrate for 5 minutes, the top silver 
will be etched away, then the nickel chrome can be etched 
away by concentrated hydrochloric acid. Rinse with de- 
ionized water, dry by dipping isopropyl alcohol, put into vacuum 
oven of 90°C overnight. 

(2) Evaporate aluminum on both sides of the wafer. 

(3) Put on photoresist: Apply photoresist on both sides by 
Waycoat 60 cps, then spin dry by 3000 RPM for 15 seconds, 
soft bake at 65°C for 20 minutes. This makes photoresist 
of 1 .4 um. 

(4) Mask making: By photographic reduction, make a mask of 

the geometry as shown in Figure 24-B. Note, that only solid- 
lined regions are for the gate-cathode mask. Some geometry 
is for the anode mask, but two stripes of 0.01 cm wide should 
be made on the edge to join these two solid-lined areas. 

(5) Expose to ultra-violet light for 10 seconds with anode mask 
and gate-cathode mask on each side of the wafers, develop 

in Kodak Thin Film Resist Developer for 30 seconds and rinse 
in Kodak Thin Film Resist Rinse for 20 seconds. 
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(6) Put into aluminum etch (one part of concentrated HC1 
and four part of deionized water) for five minutes. 

(7) Strip photoresist: Put into J-100 (at 100°C) for 10 
minutes then xylene, TCE (at 70°C), acetone for 2 minutes 
each, dry by dipping in isopropyl alcohol for 30 seconds, 
and air dry. The process is now finished. 
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APPENDIX D 


COMPUTER PROGRAMS 

Computer programs for calculating n p (x,y,t), 1 J2^ n p) * ^3* 
and Ij, are available on request from the authors. 




APPENDIX E 



I, 


*Gn 

y j 3> 





LIST OF SYMBOLS 


speed of light 

diffusion constant of electrons 
detectivity of photo-cell 
electric field 
energy band gap of silicon 
dielectric constant 
generation rate 

length of gate or cathode contacts 
anode current 

transistor reverse saturation current of the collector 
junction with emitter open-circuited 

cathode current 

gate current 

gate current of electrons 

minority carrier current in p 2 of junction J 3 

current density 

gate current density of electrons 
Boltzmann constant 
thermal conductivity 
device width 



Laplace transform of 

inverse Laplace transform of 

diffusion length of electrons/holes 

diffusion length of holes 

effective mass of electrons 

effective mass of holes 

real part of refraction index 

acceptor concentrations in p 2 -region 

effective density of states of conduction band 

effective density of states of valence band 

electron concentration in p region 

equilibrium electron concentration in p region 

= Hp (at junction J 2> when J 2 is forward biased) 

= n po ex P < t ' V J3 /kT) 

steady state concentration of n p 



hole concentration in n region 
Laplace transform of p p 
Particular solution of P 
complementary solution of P 
electron charge 



the variable in Laplace domain corresponds to time t 
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s, 2 (s + ^ )/Dp 

p 

t time 

T absolute temperature in °K 

U recombination rate 

Vj^ voltage drop in junction , i =1, 2, 3 

a optical absorption coefficient 

a.j current gain for the pnp transistor 

<*2 current gain for the npn transistor 

A wavelength of light 

u n mobility of electrons 

Up mobility of holes 

x n lifetime of electrons 

Tp lifetime of holes 

0 Theta function 

P charge density 

Pp resistivity in p region 

resistivity 

l . ; + i.; 

9x - ay.;r.., ; 
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